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ABSTRACT 

Genomic fragments of the HN and L genes from Brazilian bovine parainfluenza 3 virus (bPIV-3) isolated as 
contaminants from cell cultures and clinical specimens were amplified by reverse transcription-polymerase 
chain reaction (RT-PCR), sequenced using specific degenerate primers and analyzed by phylogenetic 
comparison with reference strains of bPI3V. The Brazilian isolates revealed a high degree of genomic when 
compared to SF4/32 prototype sfrain, within the recently proposed genotype A of bPIV-3. 
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Bovine parainfluenza virus type 3 (bPIV-3) is a member 
of Respirovirus genus of the family Paramyxoviridae, 
genetically and antigenically related to human parainfluenza 
virus type 3 (9). In cattle, bPIV-3 is associated to the clinical 
condition known as "Shipping fever", a disease characterized 
by fever, apathy, anorexia, sneezing, serous and mucous nasal 
discharges, usually related to stress induced by fransportation 
of cattle (3, 19). Occasionally, it has also been associated to 
abortions (11, 19). The virus is spread worldwide (8, 10, 20); in 
Brazil, bPIV-3 is endemically present - though seldom isolated 
- mainly linked to respiratory disorders (4, 5, 7). 

Bovine viral diarrhea virus (BVDV) is by far the most 
common and widely recognized virus contaminating cell 



cultures (1); however, other viruses may be present. Viral 
contaminants can interfere with diagnosis of viral infections, 
since its unnoticed presence can lead to interference with the 
growth of other viruses. Vaccine produced on contaminated 
cells may in turn be contaminated, leading to 
immunossupression, seroconversion or disease in the 
vaccinated animals. Other parainfluenza viruses were found as 
contaminants in cell culture (6). On the present study, bPI3V 
contamination of bovine and swine cells is reported, and the 
cell isolates were phylogenetically compared to field strains. 

In view of the limited number of available isolates, to date 
no molecular studies on Brazilian bPIV-3 were conducted. In 
the present study, the HN and L genes of two Brazilian bPIV-3 
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clinical isolates and another two bPBV cell contaminants were 
partially sequenced and submitted to phylogenetic analyses. 

Madin-Darby Bovine Kidney cells (MDBK, originally 
from ATCC, CCL-22) were used for virus production and 
assays. Cells were cultured according to usual methods (7). 
The bPVI-3 reference strain SF-4 (National Veterinary 
Laboratories Service, USA) was used as a positive confrol for 
RT-PCR tests. The bPVI-3 isolate DIO was recovered from 
nasal secretions of a calf with respiratory disease in Rio Grande 
do Sul, Brazil (7). Two other isolates were included, one from 
a case of bovine abortion occurred in the state of Goias, Brazil 
(PG1775) and another recovered from an accidentally 
contaminated culture of MDBK cells (bP13V/BR/cells). 
Another batch of Porcine kidney cells hneage (PK-15) showed 
also to be contaminated with this adventitious bPISV isolate, 
and the virus was further isolated in non-contaminated MDBK 
cells. Sfrict biosafety measures were taken on the laboratory to 
avoid cross contamination. 

All bPIV-3 isolates were cultured in virus-free MDBK 
cells following standard protocols (7). Viral RNA was 
exfracted from viral stocks with Trizol (Invitrogen, Carlsbad, 
CA, USA) and chloroform, followed by sequential 
precipitations in 40% isopropanol and 70% ethanol (15). 

Primers used were targeted to amphfy bPIV-3 genomic 
fragments as proposed by (16). For cDNA sjmthesis, viral 
RNA (1 (xg) was reverse transcribed with 200 units of M-MLV 
reverse transcriptase (Invitrogen, Carlsbad, CA, USA) and 200 
ng of sense degenerate primers. Polymerase chain reactions 
(PCR) and sequencing were performed as described elsewhere 
(18). The SF-4/32 viral strain cultured in non-contaminated 
MDBK cells was used as positive confrol and deionized water 
as negative confrol. Similarly, another PCR was conducted to 
amplify and sequence the L polymerase gene from the cell 
culture contaminants, bP13V/cells/MDBK and bP13V/cells/PK- 
15, to better infer whether they are genetically related. After 
purification and quantification, the PCR products were 
submitted to nucleotide sequence analysis. Sequencing was 



carried with the Big Dye Terminator Cycle Sequencing Ready 
Reaction (Applied Biosystems, UK) in an ABI-PRISM 3100 
Genetic Analyzer (ABI, Foster City, CA), according to the 
manufacturer's protocol. Three repetitions were performed on 
each amphfied genomic fragment. 

The sequence of SF4/32 sttain (accession number 
AF178655) and other PIV-3 sequences available at GenBank 
database (accession numbers EU277658, NC002161, 
AF178654, U31671, E01754, Y00114, E01755, AB18996, 
AB189960, AY283065, Z26523, L25350, X55803, X56131 
and X02808) were used for comparison for HN region; 
sequences AF178654, AF178655, D84095, EU277658, 
EU439429, and EU326526 were used for the L gene. The 
BioEdit software, version 7.01 (Hall 1999) was used to 
manipulate the nucleotide and amino acid of retrieved bPlV-3 
sequences. The alignment of sequences was performed using 
the ClustalW software, version 2.0 (13). 

Phylogenetic reconstructions based on the alignment of 
the HN and L genes nucleotide sequences are shown on Figure 
1. Both frees (Figure 1, A and B) were calculated using the 
Maximum Parsimony method, in a bootsfrap of 500 replicates. 
All calculations were made in MEGA5 (17). Homologous 
sequences from hPI3V were included as outgroups. 

All Brazilian cell contaminants and clinical isolates 
clearly clustered in the same group as the prototype vaccine 
strain SF4/32, within the recently proposed genotype A Since 
the 1980's, when vaccination to bPI3V was introduced in 
Brazil, there are both live and killed vaccines marketed based 
on SF4 sfrain. One may consider that the isolates presented 
here may be related to the circulation of vaccine the field, 
which may be a possibility; however, those samples were 
collected from non-vaccinated cattle or were foimd as cell 
contaminants. This, together with the nucleotide differences 
found from the original SF4/32 sttain, are indicating that 
those are wild bP13V isolates. 

Although the results presented here might be considered 
as preliminary, since they are based only in a few sequences, it 
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can be seen in Figure 1 that genotype A may present two 
distinct sub-lineages, one related to the SF4/32 strain, another 
with the 9 ION isolate; nevertheless, this have to be confirmed 
with a greater number of isolates. Analyzing the data collected 
on the present study, for the HN genomic portion, when 
isolates belonging to the same sub-lineage were compared 
intra-genotypically, a variation within the range 0.3 to 2.4% 
was recorded (with the exception of the variant isolate PG1775, 
10.6%). The overall level of inter-subgenotypical variation for 
the HN gene was higher (7.8-15%) (Tablel); however, the 
differences ranges for particular isolates do not fully support 
these conclusions, thus requiring further efforts on the 
collection and characterization of new bP13 V isolates. For the 
L gene, the identity for the analyzed genomic fragment was 
above 97% for strains inside the genotype A, and 84-89% 



between genotype B strains and viruses fi-om genotype A. It is 
also suggestive that the isolates which were noticed as 
adventitious contaminat were virtually the same virus, since the 
sequences were highly similar. Isolation of bPI3V as a 
contaminant in cell cultures may be another matter of concern 
for laboratories dealing with research and diagnosis of bovine 
viral respiratory infections, or vaccine production, since the 
cytopathic effect of the virus is sometimes very discrete. 

These findings provide additional evidence for diversity 
among bPI3V genotype A isolates, in agreement with previous 
reports based on antigenic and phenotypic analyses with 
monoclonal antibodies (2, 12-14). The impacts of this genomic 
diversity on the clinical pictures and preventive measures, 
especially vaccination, remain to be determined in the future. 
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Figure 1. Evolutionary relationship between bovine parainfluenza type 3 viruses (bPIV-3), as inferred from partial sequences of 
the gene HN (A) and L (B). The neighbor-joining method was used; only bootstrap values higher than 50 are shown. Brazilian 
strains (black lozenges) clustered together with the prototype strain SF4/32. The two distinguishable subgenotypes among bPIV-3 
genotype A isolates are shown. Human PIV-3 was used as an out-group. 
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Table 1. Identity percentages among partial nucleotide sequences (nucleotide positions 305 to 660) of the hemagglutinin- 
neuraminidase (HN) gene from different bovine parainfluenza virus type 3 (bPBV) isolates. Human PI3V was considered as an 



out-group. 
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